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1. Introduction 

Variational procedures have proved to be useful tools in the treatment of various 
quantum many-body problems, from molecular to hadron physics.'^'^] As an exam- 
ple, we shall mention the case of the low-energy hadronic spectrum, as described 
by QCE0 in its non-perturbative regime, where one has to deal, simultaneously, 
with confined fermions (quarks and anti-quarks) and bosons (gluons). To find there 
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an effective method to diagonalize the QCD Hamihonian is not at ah trivial One 
often recurs to trial states, as the coherent state defined in,*"^ which give some 
insight on the vacuum structure of QCD. IrP a first step in this direction was re- 
ported. The Coulomb interaction was approximated by a contact interaction and 
a semi-analytic solution of the QCD Hamiltonian was found. Based on these re- 
sults, we aim at developing a variational method to diagonalize the complete QCD 
HamiltoniarP within the basis constructed mW 

The paper is organized as follows: In section II a boson system is considered. 
The transformation to new quasi-particles is presented and the nature of the new 
vacuum state is discussed. In section III the same is done for a fermionic many-body 
system. In section IV we show the relation between the present method and the 
coherent state used in a variational treatment of QCD .1^1^ Finally, in section V the 
conclusions are drawn. 



2. A system of many bosons 

Let us define the creation and annihilation operators, b]^^ and b"", where the two 
indexes are to be associated with different degrees of freedom. In the case of QCD, 
the index a may refer to color and a to flavor, spin and orbital degrees of freedom. 
Lower and upper indices denote the difference under transformations of the creation 
and annihilation operators. These operators obey the commutation relation 



baa Tl 



= SalsSab ■ (1) 



We introduce the following general ansatz for the transformation to new creation 
and annihilation operators 



a. a 



V ^ p 

The indices a and (3 refer to the spatial quantum numbers, except the magnetic 
projection of total spin, and a is a short hand notation for the sum of the color 
index and the magnetic projection of the spin. By construction the matrices M 
and N transform under raising and lowering of the indices in the same way as 
the creation and annihilation operators, to preserve the transformation properties 
of the P operators. In QCD the indices a and /3 run over the values which are 
allowed by the Coulomb condition (i.e: restricted to transversal space-components) 
and they can be divided into several other indexes, as done in.l^Es' For example, 

a ^E{N,L,1)J , (3) 

where S denotes either the £ (electric) or A4 (magnetic) modes, N is the principal 
quantum number (when the harmonic oscillator is used as a basis), L the orbital 
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angular momentum and J the total spin. For the boson operators we use the phase 
convention 



(4) 



under raisirig and lowering indices, where the bar over an index denotes its conjugate 
component [^EHIll With the phase transformation ^ , we have 

Np" = (-l)v(°')+v{f» isfl^ _ . (5) 

Next, we investigate the properties of the P operators under commutation. They 
are boson operators, which gives conditions to the M and matrices. Then, 



13b 



72 fc 



j-iV M " AT h ,1 

+ 2 2.^7172 71 \" ' "72bJ 

2 \u\lia . t 

72fc 

-- , ,y ,v ''^ \b^''^°- h 

2 Z^7l72 ^'71 "'/3 " 



2 2^7172 ^^71 ^^-^0 
2 ^7l72 71 t 



'72 b 



which implies 



'5abE7l(^7"^V+^7"V) ■ 

7 

EV^7"=^"/3 • 
7 

One is free to choose Nj^ = Mj^, which we adopt from here on. 



(6) 



(7) 



2.1. The trial state for the bosonic many-body system 

For the trial state we define a new vacuum | 0), which satisfies 

P"" I 0) = . (8) 

In order to calculate matrix elements of the boson creation and annihilation oper- 
ators, Eq. ([2]) has to be inverted, giving 



(9) 



''In this notation, the angular momentum projeetion M goes into —M and in SU (3) the color- 
hypercharge and color-isospin (y, T, Tz) goes into — > (— y, T, — Tz). 
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and, from these equations one gets 

^P^-\T. {Mp"PL + M-pP^,) (10) 

The above transformation is more general. Under scaling, and preserving the com- 
mutation relations, we write 

bin V^^'ia 

6"° ^ -^6"° , (11) 

The scahng Aq may depend on the index a, a choice leading to a trial state which 
does not have definite color or spin. A further generalization of the transformation 
is possible, to include a displacement of the operators 



h' —5- \ \ h' ~ n 



The ry^a are tensors and transform in the same way as the 6-operators. 

2.2. Nature of the boson vacuum state 

With the redefinition ([Tl]) . Eq. ^ changes to 

/A^bL - -^b^a) = -V2Y,MiPpa ■ (13) 

VAa / ^ 



Applying it to the new vacuum state | 0), gives 

1 



A„6L - -4^b^a ) I 0) = . (14) 



or 

baa I 6) = A„bL I 6) . (15) 

The annihilation operators act as derivatives to the creation operators, i.e., 

d 

baa — 7— ■ (16) 
This provides, for the component (a a), the expression 

I 0) - e^^-^''"'^'"^ I 0) . (17) 



and it leads to the solution 



I 0) - e^^"^"(^'-'^'") I 0) , (18) 
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where (bl^ ■ 6^"^ = J^a^aa^""^ ■ The introduction of the parameters Aq, gives us a 
further freedom in the variational procedure. The P-operators can be cast into the 
form 



P.<. = ^>,M», |^A,s!,,8„-— !>;. I . (19) 

This allows to vary M^^ (and it complex conjugate M"^), and Xa- Using the more 
general definition of (|12l) does not change equations (|13p and (ITi|) . Thus, the struc- 
ture of the trial state will stay the same, because the differential equation in terms 
of the 6-operators will be the same. Eq. (fTO|) will change to 

A + 



V ^ a 



(20) 



and the equations for the P-operators change to 



(21) 



For the annihilation operator no contribution of rjaa enters its definition. The pa- 
rameters to vary are: 

M/ , A„ and ry„„ . (22) 

and the use of rjaa allows us to shift the operators by their vacuum expectation 
values, as done in presence of spontaneously broken symmetries. 



3. A system of many fermions 

In this section, we extend the above consideration to fermions. We use the notation 
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The b-operators refer to quarks and the d-operators to anti-quarks. They transform 
differently under SU{i). While the quarks transform with respect to color as an 
(1,0) irreducible representation (irrep), and the anti-quarks as the irrep-(0,l). The 
same transformation properties apply for the annihilation operators 

^Q(l,0)aj'm ^a(0,l)aj"m j-2^-j 

The indices a, /? refer to numbers like the principal quantum number and other 
orbital indices, except color and spin. They can be treated as cartesian. Thus, the 
phase factors (p{a) of the previous sections are set to zero, i.e. (— l)v'(") = +1. In 
the Appendix we derive the properties of the creation and annihilation operators 
under raising and lowering indices, using the phase convention of.'^'^^] Here, we 
resume the results: 

^ ' Q(l,0)a,^' — m 

_ (]\Xci+j-rn^'\a{0,l)a,j-'m 
Q(l,0)a,^'m ^ ' 

O — I,— -Lj OQ(0,l)aj-m 

b„(0,l)a,,>n = (-l)X»+^-+"b"(l'0)"'^-" , 

(25) 

with Xa = + Y + 7z , F is the hypercharge and the third component of the 
isospin. A similar property holds for the anti-particle operators, i.e., 

" — l^J-J "a(l,0)aj-m 

"a(l,0)a,jm — l^J-J " 

^to(l,0)a,jm ^ (-l)Xa+J+™dt 

\ z a(U,l}a,j— m 

dim^ - (_l)X.+.-™^t«(1.0)a.,-™ _ 

(26) 

and 

{^a(l,0)ajimi I A A A A 

" ' "/3(l,0)fcj"2m2 / ~ °''""^"J1J2"™1™2 

{^a(0,l)a,jimi »t 1 A A A A 



(all other anti-commutators vanish). To perform the mapping onto new fermion 
operators, we follow the same procedure as for These new fermion operators 
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and P are given by @ 

^a{l,0)a.,jm = X/ ^'^"^ (j^lil,0)a,jm + '^;3(l,0)ajm 

pQ(l,0)a jm ^ J_ ^.^ a 1^ ^/3(l,0)a Jm _|_ ^t/3(l,0)a 

Q(0,l)a,jm yl^ / ^ a \ p{0,l)a,jm Pl^^J-jajm 

_pQ(Oa)ajm ^ J_y^^,/a /^/3(04)a,jm ^ ^t/3(0,l)a jtt 

V2 Y 

(28) 

Lowering the indices of the annihilation operators gives (see the Appendix and 



£>Q(i,o)ajm = ^ X! ^'^a {dp(l,0)a,jni " ^^;(l,0)a,j 



(29) 



One reason for this ansatz is that there is no mixing of color (both the quark 
creation-operator and the anti-quark annihilation operator with lower indexes be- 
long to the same irrep). Note that these definitions preserve anti-commutation 
relations. 



3.1. The trial state for the fermionic many-body system 

With respect to the new vacuum, we require that 

Pa{0,l)a,jm\0) = 

^a(l,O)a,,m|0) =0 . (30) 



''For simplicity, we keep the same notation of the previous section, and the fermionic or bosonic 
character of the operators will be specified when needed. 
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By inversion we get 

^^(l,0)a,jm + <^/3(l,0)a,jm = ^ff" Pl{l,0)a,. 

a 

&/3(0,l)o,jm - dl(0,l)a,jm = '•^^ ^"^PaiO,!)^. 

a 

4(0,l)a,jm + h{0,l)a,jm = ^l3^l{0,l)a,ji 

a 



leading to the expressions 

6t -± 

"/3(l,0)a,im - ^2 
1 



Q(l,0)a,jm 



b0{O,l)a,jm - ^ E Pa{0,l)a,jm + -^1(0,1)^,, 



"/3(0,l)a,jm ^ 



D 



a(0,l)a,jrn -Pa(0,l)a,irn 



t^/3(l,0)a,jm - E -Oa(l,0)ajm + -Pa(l,0)a,jm 



(31) 



(32) 



A more general transformation can be obtained introducing a displacement, like in 
the boson case, 



a(l,0)a,jm 



^ a(l,0)a,_7m 
.J,«(ljO)a,i»n 



(33) 



for the quarks and similarly for the anti-quarks: 



it 



^" a(0,l)a,- 



a(0,l)a,jm v ck a(0,l)a,jm 

^a(0,l)a,jm ^ ^_J^^a(0,l)a,jm 



(34) 



A further generalization may be defined by introducing Grassman numbers 

^a(l,0)a,jm and r]a(0,l)a,jm SUCh 



^l(l,0)a,jTO V^^^i(l,0)a,jm Va{l,0)a,jm 

^a(0,l)a,jm ^ /r — ^a(0,l)a,jm Va(0,l)a,jm 
'A 



(35) 
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for the quark part and similarly for the anti-quark part: 

'^l(0,l)a,jm ^ y^^a(0,l)a,jm ^ Va{0,l)a,jm 

dallMaJm ^ n: — '^a{l,0)a,7"m ~ Va(l,0)a,jm (36) 



V 

The ?7Q(i_o)a.jm (^Q(o.i)aj"m) E^re tensors and transform in the same way as the d- 
and b-operators with the lower index. 



3.2. Nature of the new fermion vacuum state 

With the redefinitions (|33)) and (p4| . we obtain 



-V2E/3^4^£>/3(l,0)a,,™ . (37) 



leading to the equations 



V^'^l(0,l)ajm ~ —^ba{OA)a,jm^ |0) - 

^^i{l,0)a.jni ~ —r^d^il,0)aj7n^ I 0) = , (38) 



which imply 

^Q(0,l)ajm|6) = +-^c('^l(o,i)al^) 

da(l,0)ajm|6) = ^I(l,o)a l^) ' (^9) 

The structure of the fermion vacuum state can be determined in the same way we 
use for bosons, except for the restrictions imposed by the Pauli Principle. Denoting 
by I 0) the vacuum state of the b- and d-operators and by Af the normalization 

of the new vacuum state | 0), whose expression is IlQam ^l/-\/l + A^^ , but not 
relevant for the following discussion, we have 
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0) = AAexp 



: Afexp 



E\ jta(l,0)o,jmtt 



|0) 



= Afn^^m exp [A„dt«(^'°'"'^-'"6l(,,o)„,,„] |0) 

= n [l + (-l)'^"+'"'"^«<(0.1)a„>n^i(l,0)a. 



|0) 



(40) 



(In the last step we changed m to — m.) In the last line we expanded each expo- 
nential up to products of two terms, taking into account that j-m) ~ 

(*^a(o i)a jm) ~ This ansatz corresponds to a condensate of quark-antiquark- 
pairs, as expected. These pairs are coupled to color-spin zero, thus, the trial state 
has definite color-spin zero. Now we apply the annihilation operator 6^(0, i)a to this 
ansatz obtaining 



0) 



^ XX L a' {0,l)a'jm' a'{l,0)a'j — m' 

(a' ,a'm')^{a,am) 

= i H 1 + (-l)'^"+-' Aa'(i^,(o,i)o'jm'^l'(l,0)a'j-TO' 

y^ f_l\{Xa+j-m)+Ua+j+m)^j^a(lfi)a,j-mf^\ f_l\d\ I Q) 



n 



■^z' %'(0,l)a'jm'"a'(l,0)a'j-TO' 



^ a(0,l)ojm 



0) 



(41) 
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Since a fermion creation-operator commutes with the product of two fermion cre- 
ation operators, we have 



{a' .a' m')^ (a. am 



\ ^) "a'(0,l)a'm'"a'(l,0)a'j-m' 

xXad\^,,-. |"l + (-l)X"+^-'"Aad^„,,- . 5^,„, . 1 I 0) 

" a(0,l)ajm i ^ ' oc{0,l)a,jm a(l,0)a,j — rnj I ' 



— '^a"a(0,l)ajTO 



a'a'm' 



6) 



^ a(0,l)ajm 

We repeat the above steps with the application of the annihilation operator d 



HO) 

(42) 



(a', a' m') 7^ (a, am) 

I (a',a'm''):/(a,am) 

X (-l)2(Xa+j+r»)^^a(0,l)a,i-m^t 6^, . 10) 

^ ' " Q(0,l)a,j— m a(l,0)ajm ' ' 



I (a',a'm')7^(a,am) 



^« "a'(0,l)a'7-rn' a'(l,0)a'irn' 



X (+l)Aa&i(i,o)„,„. I 0) 



(43) 



Using again that a fermion creation operator commutes with the product of two 
fermion creation operators, we have 



da(lfl)a,3m | 0) = j\A 



n 



1 J- f — l Wo+j+m' \ J,t 

^1 '^a'"a'(0,l)a'-m'"a'(l,0)o'im' 

xAc6^ln> . [l-|-(-l)^°+^'+™Aad^„,s- . 6^, • 1 I 0) 

= +Aab|j(i o)a,jm j II \^ ^ ^'^a' {0,l)a' jrn'^^a' (lfi)a' j-ra' \ 

\ <x'a'm' ) 

= -VXJb^ . 10) . 

This proves that our ansatz satisfies the operator equation for the new fermionic 
vacuum. The vacuum has definite color and spin, when we assume that the indices 



|0) 
(44) 
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of the transformation matrix do not depend on color nor on the spin quantum 
numbers. With this, the P and 7D-operators are given by 



0)a,jm 



Adding a displacement (by adding a Grassman variable), the equations for the 
P-operators change to 





Pa{0,l)ajm = ^''a ^ y=^>/3(0,l)a,jm " \/^d^fHo,l)a. 

,l)a,jm 



(46) 



In total, we have now the following parameters to vary: 

Mj , Xa and Va{l,0)a,jm (47) 

Finally, wc give the expressions of the former creation and annihilation operators in 
terms of the new ones, including the information on Aq,, r]a(i,o)a,jm and ??a(o,i)a,jm- 
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P(l,Q)a,jm 



5/3(0, l)o,im 



1 



■13 



A/3??/3(0,l)a,: 



1 



'7/3(0, l)aj"m 



-^1(0, 1)3, jm -fQ(0,l)a,im 



/A 



•/3 



d;3(l,0)a,im " V E ^'^Z?" -Da(l,0)a,im + -Pl(1^0)a, 



(48) 



4. Relation to the coherent state as used irJ^I^ 

IrPI^ a monopole gas for the description of the QCD ground state was proposed 
and thermodynamical properties extracted. The trial state of this monopole gas, 
depending on several parameters, has the form of a coherent state. In this section we 
show that this trial state can be recast in the many-body language of the previous 
sections. The trial state introduced irPEl has the form 



where A/" is a normalization. We start with the expression of the field components 
A{x)a in the basis of the operators b' and b 

A(x)„-^0"(x)^(bL + baa) 

At (x), = A'^ [x] =J2K (x) ^ {b^'^" + b"-^) , (50) 
The amplitude (p^{x) is to be associated with the factor of,^ and 

<P*Jx)^r{x) = (-l)^("Va . (51) 

In the quantization of the field Aa{x) there appears an additional factor l/^/TlZ, 
with being the frequency of the solution (j)°'{x). We adopt the notation that this 
factor is included in the function. In order to preserve gauge invariance,'^ (x) has 
to contain a combination of perturbative A-fields with a monopole solution, i.e., in 
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general, different boson operators have to be defined. For simplicity, we use one 
generic expression. 
With this, we have 

Al{x)uj{x - y)Aa{y) = 

71 Ea^ <j)aix)uj{x - y)(j)'^ {y) 

(6"'^6t^ + 6t-6,, + 6t-6t^ + 6"''6,,) . (52) 
where, in lowering indexes and contracting them we have used the relation 

= r{x)b% , (53) 
With this, ([52]) can be rewritten as 

-j^EaJ:^pr{xMx-y)cf>^{y) 

{blb\^ + btbp, + btbl^ + blbf^a) ■ (54) 

In terms of the scalar products, 

Y.blb\^^{b^-bl) . (55) 

a 

it reads 

Y.^Al{xMx-y)Aa{y) = 
7lEa/^r(a:)w(a;-y)0'3(y) 

{ (ba • 6t ) + (fot . + (bi-bl)+{b^-bp)] . (56) 

Because a is a dummy index, we can skip the bar over a in the sum. 

By integrating Eq. (j52p over x and y we obtain the same structure as in (1561) . 
by writing 



^"'3= I dxdy(l)°'{x)uj{x-y)(j)'\y) 



(57) 
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for the factors in front of the scalar products. Using the ansatz for monopoles given 
inPI^the factors A"-^ can be determined. What we have shown is that the coherent 
trial state, as used in,'^'^ can be cast into a standard language, like in conventional 
many body theories. 



4.1. Taking the product of the operators P' and P 

The expression (j54p for the exponent can finally be rewritten as a product of cre- 
ation and annihilation operators by taking the product of two P-operators as 

[iVI^b^^" - N'^pb'^^^"-^ . (58) 
Performing the explicit multiplication gives 

+Bf;{b,,-bP^)-Af:{b,,-b^^^)} , (59) 

with 



a 
a 
a 

^h=\T.^o.Nj^n. ■ (60) 

a 

In order to get the same structure as in ((56)) - ([57]) . the upper indices in the b 
operators in (l60l) are lowered and the corresponding lower indices in the coefficient 
matrices are raised. In addition, the canonical transformation bj^^ — >■ ibj^^, — >■ 
is performed, which maintains the commutation relations. This transforma- 
tion changes the sign in front of the products of two creation or two annihilation 
operators, but not in the product of a creation with an annihilation operator. With 
this, we have 

+B'P^P^ {bp, ■ -I- A"^^l'^ (bp, ■ bt J } , (61) 
We choose the particular relations 
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thus, leading to ([56| with the integration ([57l) performed. ((62|) demonstrates that the 
ansatz of a monopole gas is contained in the many-body trial state, but eliminating 
the conditions ()62p permits a more general structure, which is of great advantage. 
The QCD trial state can now be formally written as 

Because P"" | 0) = 0, only the first term in the exponential contributes. Thus, the 
trial state has been reduced to the vacuum state \ 0), 

The matrix elements MJ^ (= NJ^) can be dealt with in two ways: 

i) Simply as parameters which are determined by minimizing the expectation value 
of the Hamiltonian with respect to a trial state, or 

ii) determine Mj^ via the A"^ of Eq. ([57)1 using also the relations (15(11) - (| 6 2 1) , i.e., the 
matrix elements depend then on the parameters of the monopole gas. The above 
relation may be affected by spin mixing, because when the indices a and (5 contain 
the total spin j, the matrices will mix the spin and the trial state does not 
have definite spin, though the new vacuum may have definite spin. If we want to 
avoid this, a and /3 should not depend on the total spin. 

5. Conclusions 

In this contribution we have developed a variational procedure which includes si- 
multaneously fermions and bosons. We use explicitly the case of non-perturbative 
QCDjI^f^as a working example, but it can be generalized to any system of fermions 
and bosons. A series of trial states have been proposed, with increasing complexity. 
They include a simple unitary transformation plus a possible re-scaling of the boson 
and fermion operators and a shift of the boson operators. 

We have shown that the trial state presented in,'S]El where a coherent state 
written in terms of a QCD functional was used in order to describe the ground state 
of QCD as a monopole gas, can be recast into a standard many-body framework. 
We think that this connection may facilitate the use of such techniques in dealing 
with the low energy domain of QCD. Further work is in progress concerning the 
use of effective QCD-inspired Hamiltonians. 
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Appendix 

In this appendix we derive properties of the fermion creation operators under low- 
ering and raising their indices. We start from the convention we have introduced 
inP which is 



j^t-«a(0,l)a,jm ^ (_ 1 ) 5 +? (_ 1 )Xa+J -m fet 

^ / ^ / -\-^a{l,Q)a.j—m 

(64) 

(a, /3 correspond to Cartesian indices). In the above equation (|64|) . a is the color 
index, j is the spin, m its projection and ^ is the pseudo-spin component. The 
pseudo-spin quantum number refers to the upper and lower level and for the upper 
level ^ = 5 and for the lower level it is ^ = — ^. The phase factor Xa is always 
integer, while j and m are half integer. 

The property of the annihilation operator, under raising and lowering indices, 
is obtained by the hermitian conjugation of ()64p 

/'5t-«a(0,l)ajmy ^ ^ 



-^a{0.1)a,jm 
^ ' ^a(l,U}a,j— m 

(65) 

Identifying the fermion creation operator with component ±i with the creation 
of a particle (b^) and the annihilation of an anti-particle ( d), respectively, and sim- 
ilarly for the fermion annihilation operator , we arrive at the following expressions 



^ ' a(l,0)a,j— m 

a{l,0)a,jni V / 

" — { ^) "a{QS)a,j-m 

h , , — ( ^\Xa+3+mla{l,0)a,j-m 

(66) 

A similar property holds for the anti-particle operators, i.e., 



^a(0,l)o,im _ ('-nXa+J-mj 

" — V "'a(l,0)aj— m 

"a(l,0)aj'm — l^-^j " 
jta(l,0)a,jm _ 1 ^Xa+i+n^ jt 

dliO.iya.rn = (- l)>^«+^-™dt"(1^0)-^-™ 

(67) 
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Here, one has to include an additional change in sign due to the phase (—1)^^^. 

Note, that the resuhing phase property is the same for the particle creation as for 
the anti-particle creation operator. The same holds for the annihilation operators. 

With the help of the just obtained results, we can study the structure of the 
D- and D'^-operators (see Section III). The creation operator corresponds to 
a particle-annihilation operator at negative energy. Thus we make a similar ansatz 
for it, as for the annihilation operator f , but with a different sign, i.e., 

I /3(l,0)a j' 

" ^2 Y ' ^ ^ 

(68) 

It can be easily shown that this anti-commutes with f ^. Now, we lower the index 



in ([68)) . leading to 

^ (_l)Xa+j+mjr)t 



j-m 



^i3{a.l)a,j-m 



(69) 



The matrices M satisfy M/ = M^. From this we conclude that the correct ansatz 
for the new anti-particle operator is 



A similar manipulation for the anti-particle annihilation operator 

Da{l.,0)a,jni = ^ E ^^"^ (c^/3(1.0)a jm - ^')3(l,0)aj: 



(70) 



yields 



j-fa{0,l)a,jm ^ J_ f ^/3(04)ajm , ^t/3(0,l)a jm\ 



(71) 



(72) 



and 



£'a(l,0)a,,™ - (_l)X.+J+™I3"(0.1)a.J-™ . (73) 
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Equivalent considerations are applied for the raising and lowering of the indices 
of the and P operators, with the result: 
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